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Abstract 



We present new one loop calculations that confirm the theorems of Joglekar 
and Lee on the renormalization of composite operators. We do this by con- 
sidering physical matrix elements with the operators inserted at non-zero mo- 
mentum. The resulting IR singularities are regulated dimensionally. We show 
that the physical matrix element of the BRST exact gauge variant operator 
which appears in the energy-momentum tensor is zero. We then show that 
the physical matrix elements of the classical energy- momentum tensor and the 
gauge invariant twist two gluon operator are independent of the gauge fixing 
parameter. A Sudakov factor appears in the latter cases. The universality of 
this factor and the UV finiteness of the energy-momentum tensor provide an- 
other method of finding the anomalous dimension of the gluon operator. We 
conjecture that this method applies to higher loops and takes full advantage 
of the triangularity of the mixing matrix. 



PACS number(s): 11.15.Bt,12.38.Bx 
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I. INTRODUCTION 



The anomalous dimension of the twist two gluon operator M is an important quantity 
in QCD phenomenology and has been the subject of study for many years. It contributes 
to the logarithmic violation of Bjorken scaling which was one of the great early successes 
of QCD. Its inverse Mellin transform is the Gribov-Lipatov-Altarelli-Parisi gluon splitting 
kernel. In the late 1970's two groups independently calculated the anomalous dimension of 
the gluon operator in second order. One group worked in the axial gauge and the other in 
the Feynman gauge ||. The results were different, and the questionable axial gauge result 
was chosen to be correct on the basis of a relation among anomalous dimensions derived 
from supersymmetry [|J. The third order anomalous dimension of the gluon operator is 
presently needed to complete the second order QCD corrections to the deep inelastic struc- 
ture functions F 2 (x, Q 2 ) and Fl(x, Q 2 ) |J. In preparation for this third order calculation one 
must completely understand the discrepancies between results of the two loop calculations. 
Therefore Hamberg and van Neerven |J recently recalculated the anomalous dimension of 
the gluon operator to second order in the Feynman gauge. Their result agreed with the 
previous axial gauge result. However, the situation is still not entirely satisfactory, as it is 
not obvious that the new calculation satisfies the general theorems on the renormalization 
of composite operators in gauge theories J7| , ■ 

There are three points of contention. First, the physical matrix elements of the gauge 
variant (GV) operators || used by Hamberg and van Neerven do not vanish. Second, a 
physical matrix element of the gauge invariant (GI) gluon operator turns out to be gauge 
dependent. Third, their renormalization mixing matrix is not triangular. While these find- 
ings sound bad, they do not contradict the general theorems J7|, which state that the 
operators should be inserted at non-zero momentum. Since Hamberg and van Neerven 
worked with operator insertions at zero momentum the general theorems do not apply. This 
point has recently been emphasized by Collins and Scalise [|l(J who isolated an infra-red (IR) 
singularity in the proof that the physical matrix elements of BRST exact operators vanish. 

To clarify the situation we consider it necessary to verify the general theorems by explicit 
calculations. To our knowledge this has never been done mostly because there was no clear 
way to handle the IR singularities encountered when putting the external legs on-mass-shell 
0, [Q. The on-mass-shell renormalization scheme is not well defined in QCD but one may 
try a variant of this scheme wherein the external legs are put on-mass-shell in the limiting 



sense only. This leads to the so called "modified LSZ" prescription ||10|| . We show below 
that if one tries to verify the gauge parameter independence of a physical matrix element 
of the gauge invariant classical energy-momentum tensor inserted at non-zero momentum 
using this method the gauge dependence remains and the result contains logarithms of the 
external momentum squared. Hence the IR singularities, which arise when putting the 
external momenta on-mass-shell, spoil the gauge independence of the result. We present a 
calculation in which we resolve this problem. We have found, by inserting the GI operators 
at non-zero momentum and using dimensional regularization to regulate the IR divergences, 
that there are universal double and single IR pole structures in the physical matrix elements 
of both the energy momentum tensor and the gluon operator. Performing the calculation 
in this way in space-time dimension n > 4, we verify the general theorems. By appealing to 
the ultra-violet (UV) finiteness of the energy-momentum tensor, we find that the difference 
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in the coefficients of the single pole terms is precisely the anomalous dimension of the twist 
two gluon operator. This approach is very similar to the factorization of mass singularities 
in perturbative QCD and the way one deals with Sudakov logarithms. A review of such 



logarithmic factors is given in [[LI]]. We believe this method will also work for the anomalous 
dimension in higher loops and an all-orders proof of the IR factorization should be possible. 

This paper is organized as follows. In Section II we introduce our notation for the 
operators. The Feynman rules are given in Appendix A. In Section III we list some of 
the predictions of the general theorems on the renormalization of composite operators. In 
the main section, Section IV, we present an explicit first order calculation of the physical 
matrix elements with the energy-momentum tensor and the twist two gluon operator inserted 
at non-zero momentum using both of the methods introduced above. Appendix B contains 
information on the evaluation of the integrals needed to perform the calculations. We discuss 
our findings and give our conclusions in Section V. 



II. NOTATION 

We consider the effective Lagrangian 

C = Cym + C g f + Cfp , (2.1) 

where Cym, C g f, and Cpp are the Yang-Mills, gauge fixing, and Faddeev- Popov ghost 
Lagrangians, respectively, and are given by 

Cym = -\f;„F^\ (2.2) 

In this paper we denote renormalized quantities by the subscript r. For example, A^ r is 
the renormalized Yang-Mills (YM) field and A ^ is the bare YM field, u a is the bare anti- 
commuting ghost field, and £ a is the bare anti-commuting anti-ghost field. The non-Abelian 
field strength F£ u and the covariant derivative D ab are defined in terms of bare quantities as 

F« v = dpAZ - d v Al + gr hc A\Al , (2.3) 
Df = 5 ab d,-gf abc A^. 



The structure constants of SU(3) are defined by [T a ,T b ] = if abc T c where T a are the gen- 
erators of 577(3). The bare gauge fixing parameter is a, and the bare coupling constant 
is g. The YM Lagrangian is (constructed to be) invariant under the infinitesimal gauge 
transformation 

6A a ^ = --Dfe b , (2.4) 

g 



while the effective Lagrangian (2.1) is invariant under the infinitesimal BRST | 1~2| , fl"3 
transformation 
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8 BRST Al = \Dfu\ (2.5) 
with A anti-commuting, provided the ghost and anti-ghost transform as 

5 BRST uj a = - l -gf abc uj b u c \, (2.6) 
a p 

The energy-momentum tensor || follows from the effective Lagrangian by differenti- 
ation with respect to the metric tensor |T5[] without use of the equations of motion. It is then 
contracted with two powers of an arbitrary, light-like vector A M to pick out the symmetric 
traceless part. We define 

o em = -\^^e, v = oZ + oZ , (2.7) 

where 

OZ = -~A<W0gf = iA^A'F-' , (2.8) 
OZ = ~\ A^JT , (2.9) 

= --A ■ A a A ■ d{d ■ A a ) - (A • <9£ a )A • D ab u b . 
a 

The Feynman rules for OZ an d OZ which we require are given in Appendix A. One can 
show that OZ is invariant under the gauge transformation defined above and that OZ is 
BRST exact, meaning that it can be written as a BRST variation of another operator called 
its ancestor. Explicitly, 

oZ = s B R ST n, (2.io) 



with 



n= (A-3f)A- A a , (2.11) 



as can be checked using (2.5) and (2.6). 

In physical applications of perturbative QCD one is interested in the anomalous dimen- 
sion of the twist two gauge invariant gluon operator [l|]. At zero momentum, contracted 
with m powers of an arbitrary light like vector, A^, it is 

O^(0) = —P ai [) aia2 . . . £) a m-2a m -i pa m -ia ^2 12) 

g 2 a 

where 

pa — \fj, ™ 
a a/i ' 

D ab = iA H D ab 

(J, 
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and m = 2,4,6, . . . . One can show that Og is invariant under the infinitesimal gauge 
transformation defined in (2.4). We write Og schematically as F(D m ~ 2 F). As we intend 
to insert GI at non-zero momentum the operator to be considered is 

Of{Q) = F(D m ~ 2 F) + (DF)(D m - 3 F) + ■■■ + (D m ~ z F)(DF) + (D m ~ 2 F)F . (2.13) 

One can show by explicit construction that it is possible to find a A M such that A • Q = 
but Q = pi — p 2 7^ 0. So that 

S gi = J d^xe ix - Q G g \Q) , (2.14) 

can be integrated by parts and put in the form of the zero momentum case above, i.e. 
F(D m ~ 2 F). For practical calculations this involves nothing more that using the relation 
A • pi = A • p2- For convenience, the normalizations of and Og are chosen so that 
at m = 2 they coincide. The Feynman rules for GI which we require are also given in 
Appendix A. 



III. GENERAL RESULTS 

A large body of literature exists on the renormalization of composite operators in gauge 
theories. Various gauges have been explored including the axial, covariant and background 
field gauges. In this paper we will only be concerned with the covariant gauge results. 

The covariant gauge investigations began with the work of Gross and Wilczek [|1|] and 
culminated in a series of papers by Joglekar and Lee, and Joglekar 0. Most intermediate 
references are cited in |?J with the exception of || where the first order UV pole pieces 
of the energy-momentum tensor were calculated in the general covariant gauge and shown 
to add to zero after wave function renormalization. This calculation was done for non- 
zero momentum insertion and showed that the energy-momentum tensor for pure unbroken 
Yang-Mills is finite to first order. However, no effort was made to retain the finite pieces 
(to see if they were gauge independent) nor were the physical matrix elements of the Of^ 
operator calculated (to see if they vanish). Also of interest is the simplified proof of some 
parts of the Joglekar and Lee paper using cohomology techniques |L6[ and a recent detailed 



study of the energy-momentum tensor inserted at zero momentum |L0| wherein Collins and 
Scalise isolate an IR divergence that spoils the proof of Theorem 2 below. 

We begin with some terminology. A non-zero momentum insertion means that there is 
a momentum flowing through the operator vertex. For example, we take pi flowing out one 
leg, P2 flowing in the other, and Q = p\—p2 flowing into the vertex as shown in Figure (la). 
Gauge independent means that the result is independent of the gauge fixing parameter a r . 
Gauge invariant means that the operator is invariant under the gauge variation defined in 
Eq. (2.4). A physical matrix element of an operator is an on-shell truncated proper (on- 
shell amputated one-particle-irreducible) Green function with an insertion of an operator. 
By on-shell we mean external legs are contracted with physical polarization vectors, put 
on-mass- shell, and properly renormalized. 



We now quote without proof some of the theorems M, M, [17] on the renormalization 



of composite operators in non-Abelian gauge field theories. As has been noted in 
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and [1(J these theorems only apply at non-exceptional momentum. That is, the operators 
should be inserted at non-zero momentum. 

Let Gi be a set of gauge invariant operators which mix among themselves under renormal- 
ization, Aj be the set of gauge variant operators which mix with Gi under renormalization, 
and Ei be the set of operators which vanish by use of the equations of motion which mix 
with d under renormalization. 

Theorem 1 A» can be chosen so that it is BRST exact up to equations of motion, i.e. 

Ai w SBRsrBi . (3.1) 

Theorem 2 Physical matrix elements of BRST exact operators are zero. 
Theorem 3 Physical matrix elements of Ei are zero. 

Theorem 4 Physical matrix elements of Gi are independent of the gauge fixing parameter 

Theorem 5 The renormalization mixing matrix is triangular, i.e. 

(R[G}\ ( Z GG Z GA Z GE \ [G\ 

R[A] = Z AA Z AE \ A . (3.2) 
V R [E] ) V Z EE ) \E ) 

It is Theorems 2, 4 and 5, that, at first sight, appear to be violated by the work of Hamberg 
and van Neerven |J who, following [|IJ and ||, worked with zero momentum insertions. We 
demonstrate in the next section that these theorems hold to one loop when one considers 
non-zero momentum insertions and works in space-time dimension n > 4. 

IV. EXPLICIT CALCULATIONS 

In this section we outline the calculation of the first order physical matrix element with 
two external gauge fields and an insertion of Of^, Of^, and O^ 1 at non-zero momentum. 
We do this in two ways. First, we try a traditional LSZ prescription and a "modified LSZ" 



prescription [10] and see that these methods fail to respect Theorem 4 given above. Second, 
we work in n > 4 space-time dimensions and consider truly massless external legs and find 
that the theorems given in the previous section are confirmed in our one loop example. As 
mentioned above, the first order pole pieces of the matrix elements with an insertion of the 



energy-momentum tensor at non-zero momentum have been studied in detail in ||. In [10 
a detailed study of the energy-momentum tensor inserted at zero momentum was made. 
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We calculate matrix elements of operators inserted at Q = p\ —p2 ^ thereby introducing 
another mass scale Q 2 into the problem. The algebra is therefore more involved than in the 
case of Q = 0. The tensor integrals were reduced to scalar integrals as discussed in Appendix 
B. For this process we performed all algebra on a computer using the program form ||18|| . 



except for the inversion of the large matrices involved, for which we used maple v |19| . The 
Feynman rules for the vertices containing 0^, O^, and O^ 1 are given in Appendix A. We 
used the gauge field propagator 

^(*) = ^(^-(l-«r)^) > 

in all the calculations which follow along with the other standard vertices from the effective 
Lagrangian. 



A. Calculating in n = 4 dimensions 



In calculating physical matrix elements when truncating the external propagators one 
keeps a square root of the residue of the pole of the full propagator for each leg as one goes 
on-shell. The other square root is discarded with the truncated propagator. To find the 



residue, note that the full propagator is p0| 



6 ab 
I 2 



Qnv k ll k u /k 2 kak v 
+ a, ■ 



i + n(fc 2 ) 



k 2 



(4.2) 



The self energy of the YM field is well known [2TJ. It is derived by expanding terms (—k 2 ) 6 ^ 2 
for small e, using x e = 1 + eln(x) + • ■ ■, then after the contribution from the Lagrangian 
counterterm is added the 1/e terms cancel and e is set to zero. However, k 2 is not set to 
zero due to the presence of terms in ln(— k 2 ). We write 



Uf u (k) = 5 a \k^k u - k 2 g^)Tl(k 2 



(4.3) 



then 



U(k 2 ) 



9r 
167T 



1 -k 2 
l + iln(^) 



a. 



(4.4) 



where we have used the MS scheme. The wave function renormalization constant in this 
scheme is 



97- 



In the above equations e 



7-\ 



.-1 



16tt'- 



;C A - 



10 



+ (1 



iln47r, e 



(4.5) 



n 



4, f acd f bcd = 5 ab C A , 7 is the Euler 



co6 



constant, and we have extracted the dimension of the coupling by introducing the mass scale 
fi. We are now in the position to extract the coefficient of the pole as k 2 — > 0. It is normally 
chosen as 
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k+p, 



1 
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\ 
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5? 



1% 
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Q=Pr 



(a) (b) (c) (d) (e) (f) 

FIG. 1. Diagrams contributing to the physical matrix element of OfJ t and O^ 1 to one loop. A 
circle with cross denotes an operator insertion. 



R 



l + n(Ar%^ ' 



-i 



1 - 



_9r_ 

167T : 



_ f 31 5, -k 2 



+ (1 



1, -k\ 



4(1 -«r) 



(4.6) 



We will see below that if we use 



g[Pn(fc 2 )] 

dk 2 



-1 



# r J 31 + 5 n 
T6V CA rT + 3 (1 



ln(-^)) + C 1 



av 



1. 

2 /i 2 



i' 1 



av 



(4.7) 



as the residue, as suggested by Collins and Scalise [ItJ, then at Q = the physical matrix 
element of O em = + Of^ remains unaffected by higher order corrections and is indepen- 
dent of the gauge fixing parameter a r . This absorbs infrared singularities into the asymptotic 
in and out states much like what one does in proofs of factorization in perturbative QCD. 



1. Energy-Momentum Tensor 



a. Gauge invariant operator, O^. We begin by considering the gauge invariant piece 
of the energy-momentum tensor. The diagrams relevant for the calculation of the two 
leg physical matrix element at one loop are shown in Figure 1. The one loop amplitude 
A ab (— pi, p 2 , 2) is given by diagrams (b) - (e). Diagram (f), which is constructed using the 
four leg 0^ vertex is zero in dimensional regularization so it is not considered. We were 
careful to include the symmetry factor 1/2 for diagrams (c), (d), and (e) and then partial 
fraction the resulting expression to remove all terms of the form k ■ p\ and k ■ p2- Next we 
contracted with polarization vectors and used the on-shell conditions e(pi) -p\ = = e{j>2)-P2 



and p\ 



m 



but refrained from setting m = 0. 



As a check, we set Q = 0, performed a reduction to scalar integrals, and expanded the 
resulting expression in powers of e = n — 4. The physical matrix element is, 



< O%{0) > = RZ 3 



T, 



ab i 



p,p,2)+Af v (-p,p,2) e»( P y(p) 



(4.8) 



-6 ab (A-p) 2 + 6 ab -^C A (A-p) 2 

l07T Z 



3 1 

-(1 - Or) + ^(1 - Or) 
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(a) 



(b) 



(c) 



(d) 



® 

/ \ 



KO J,' 



/ V? 

/ KO 
I I \& 

J I \o 

* I \£> 

1 j 

1 7 \o 



i w 



(e) 



(f) 



(g) 



FIG. 2. Diagrams contributing to the physical matrix element of Ogm to one loop. Momenta 
routings are the same as Figure 1. 



where we have used pi = p 2 = p, e(p) • p = 0, p 2 = 0, and e(p) ■ e(p) = — 1. If we would have 
used R instead of R the coefficients of {1, (1 — a r ), (1 — a r ) 2 } would have been {4/3, —2, 1/4}. 
This agrees with || and QTUJ and provides a check on our calculation to this point. 

Having made this check we now proceed with the Q 7^ calculation. The only difference 
is that the algebra and integrals become more involved. A brief discussion of the methods 
is found in Appendix B. We quote our findings for the physical matrix element of as 
follows: 



< M2 (-Pi,P2,2)l e^(p 1 )e^(p 2 ) 



< 0^(Q) > = RZ 3 [T^ 2 (- Pl ,p 2 ,2 

= e^( Pl )e^(p 2 )T^ 2 (- Pl ,p 2 ,2) 
- 2 r ll 11, (Q 2 



(4.9) 




3 (A-p 1 ) i -(A-p 1 )(A-p 2 ) + -(A-p 2 )^ 

It is clear that this expression is ill defined at m 2 = (so there is a problem with infrared 
divergences) and is dependent on the gauge fixing parameter a r . 

b. Gauge variant operator, Of^. We now consider the gauge variant piece of the energy- 
momentum tensor. The diagrams relevant for the calculation of the two leg physical matrix 
element at one loop are shown in Figure 2. The momenta are routed similar to those in 
Figure 1. The tree level 0^% diagram (a), does not contribute as it vanishes when contracted 
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with polarization vectors. The one loop amplitude B^ u (— pi,p2) is the sum of diagrams (b)- 
(g) including the symmetry factor 1/2 for diagram (c). We partial fractioned the resulting 
expression and used the on-shell conditions as was done for the GI case. 

At this point, as a check, we set Q = 0, performed a reduction to scalar integrals, and 
expanded the resulting expression in powers of e = n — 4. The physical matrix element is, 



<O^(0)> =RZ 3 V£(-p,p) + B»(-p,p) e"(p)e"(p) 



(4.10) 



B 



ab i 

fJLV \ 



l07T z 



3 — -(1 — a r 



i (l-«,) 2 



where we have used p\ 
with [10 1 and gives 



p 2 = p, e(p) ■ p = 0, p 2 = and e(p) • e(p) = —1. This again agrees 



< O em (0) >=< o£(0) > + < 0^(0) >= -f(A.p) 2 



GV, 



(4.11) 



The form of the physical matrix element of O em is unaffected by higher order corrections 
and independent of the gauge fixing parameter with the use of R as the residue instead of 
R. 

We now proceed with the Q non-zero result. The only difference is that the algebra and 
integrals become more involved. For the physical matrix element with inserted and Q 
non-zero we find zero. Explicitly, 



< oZ{Q) > 



RZ 3 



V, 



ab 



WM2 v Pl,P2)+B- W2 {- Pl ,p 2j 

-Pi,p 2 )e fll (pi)e^(p 2 ) , 



ab 



'2) , 



(4.12) 



0. 



Which is in agreement with the Theorem 2 and is not affected by what we choose for R 
because the tree level diagram V^ l „ 2 (— P\,P2) vanishes when contracted with polarization 
vectors satisfying e(p) • p = 0. 

By looking at Eq. ( f4.8p we see that Theorem 4 is violated so at this point we already 
notice a problem. By looking at Eq. (|4.10|) we see that Theorem 2 is also violated. However, 



Eq. ([4.12 ) confirms Theorem 2 but Eq. ( [4.9|) has IR divergences in the form ln t (Q 2 /m 2 ), 



They are due to the order of our limits, namely, we took the 



1, 2 that must be removed, 
limit e — >• too soon. 

Since there is an IR divergence we can extract it by reworking the calculation and taking 
the limits p\ — > and p\ — > before e — > 0. In dimensional regularization, with space- 
time dimension n = 4 + e, both the UV and IR divergences are regulated. After UV 
renormalization, which is performed at e < we can analytically continue to e > which 
regulates the IR singularities. This allows us to set the square of the momenta of the external 
gluon lines to zero because (pfY^ 2 = = (p 2 ) e ^ 2 so long as e > 0. This interchange of the 
order of the limits is the crux of our new calculation which we present next. 
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B. Calculating in n > 4 dimensions 

In contrast to the previous section, here we reverse the order of limits and find that 
Theorem 4 is now satisfied. That is, we consider e > so that (— k 2 ) e ^ 2 = in the limit 
k 2 = 0. This implies that the residue factors RZ 3 of the previous section are to be set to 
unity. We will see that doing this leads to a universal double and single pole structure that 
allows use to determine the anomalous dimension of the gluon operator. 



1. Energy -momentum tensor revisited 



a. Gauge invariant operator, 0^ m - We now return to the gauge invariant piece of 
the energy-momentum tensor. The diagrams relevant for the calculation of the two particle 
physical matrix element at one loop are shown in Figure 1. Diagram (f), which is constructed 
using the four leg vertex, is zero in dimensional regularization. Also, because we set 
(p 2 ) e ' 2 = = (plY^ 2 diagrams (d) and (e) vanish. We partial fractioned the expression 
resulting from diagrams (b) and (c) to remove all terms of the form k ■ p\ and k ■ p 2 . Next 
we contracted with polarization vectors e w (pi) and e^ 2 {p 2 ) and used the on-shell conditions 
e (Pi) " Pi — — e G°2) ■ P2 and Pi = 0=P2 remembering that (p\Y^ 2 = = (p 2 Y^ 2 as e > 0. 
A discussion of the methods for evaluating the resulting integrals is found in Appendix B. 
Adding diagram (a) we quote our findings for the sum of all the diagrams in Figure 1 with 
inserted at Q and external legs on-shell as defined above as follows: 



<0^{Q)> =1 T;l 2 (- Pl ,p 2 ,2)+A^ 2 (- Pl ,p 2 ,2) e^( Pl )e«(p 2 ) 



e^( Pl )e^(p 2 )T^(- Pl ,p 2 ,2) 



x il+g' r C A S n (- 

+ glCASni — ^J 



ab 

Q 2 . 



8 -4 7g + 22/3 + /i(2) 



e(Pi) ■ e (P2) - 



Pi ■ e(p 2 )p 2 ■ e(pi) 
Pi -P2 



x 



l(A.p 1 ) 2 -(A.p 1 )(A.p 2 ) + ^(A.p 2 ) 2 



(4.13) 



where Ca = N c = 3, S n 



71 



n/2 



/(2tt)' 



and T^ 2 (- Pl ,p 2l 2) 



is the Born level vertex for 



the classical energy-momentum tensor given explicitly in Appendix A. The factor of pj 



-2c 



is the dimension of the coupling g r . The finite term h(2) is given later in Eq. ( 4.17 ) for 
completeness. 

At this point we pause to note several important things. One, this result is independent of 
the gauge fixing parameter in agreement with Theorem 4. Two, because we are considering 
on-shell legs we have not multiplied by an LSZ type residue when truncating the legs. Three, 
there are both double and single Sudakov IR poles. Four, there is a new finite form factor 
that was not present at the Born level. Five, if we take the limit Q 2 — > for e > the only 
piece that survives is the Born level vertex. We will return to these points in the discussion 
below. 
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b. Gauge variant operator, Of^ . We now return to the gauge variant piece of the 
energy-momentum tensor. The diagrams relevant for the calculation of the corresponding 
physical matrix element with an insertion of at one loop order are shown in Figure 
2. The momenta are routed similar to those in Figure 1. Again, (f) and (g) vanish. The 
tree level Of^ diagram (a), does not contribute as it vanishes on-shell. After performing 
the integrations as above the result vanishes, which is in agreement with the Theorem 2. 
Explicitly, 



< o^Z(Q) > = i [v^ 2 (- Pl ,P2) + b;1 1 (- PuP2 ) e^( Pl y>( P2 ) , (4.w) 

r>ab , 

= 0. 



B^(- Pl ,p 2 )e^( Pl )e^( P2 ) 



2. Gluon operator 

Now consider the twist two gauge invariant gluon operator. The one loop diagrams are 
shown in Figure 1. We followed the same procedure as outlined above for the revisited 
energy-momentum tensor but this time we took the liberty of setting A • p\ = A ■ p 2 = A ■ p 
which allows us to do the integration by parts described at the end of Section II. The 
necessary integrals are discussed in Appendix B. We quote the result for the sum of all 
diagrams on-shell as defined above: 



< Of'{Q) > = 1 



T^ 2 (- Pl , P2 ,m) + A a J> ifM2 (-p 1 ,p 2 ,m)\ e^( Pl )e^(p 2 ) 
(pt)e^(p 2 )T^ 2 (- Pl ,p 2 ,m) 



(4.15) 



x h+£c A S n (-^y -- + 
+ g ^C A S n (-%y(A-p) m f(m 



8 -4 7B + 22/3 + 7 J?VCU 



+ h(m) 



e(Pi) ■ e fe) - 



P i ■ e( P2 ) P2 ■ e(pi) 



P i - P 2 



where T^^ 2 (— P i, P2 , m) is the Born level vertex for the gauge invariant gluon operator given 
in Appendix A. All other symbols are as defined below Eq. (|4.13|) and 7^ is the one loop 
anomalous dimension of the gluon operator H 
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+ 



16 



771—1 



m 



rn 



2 m + 1 

The finite functions h(m) and f(m) are given by 



m — 



1 i=i 



h(m) 



-7l + 



jEhJt ,22; 
2 I c A ^ 3 , 
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8S(m-2) -SS^m 
8 12 



I! 



22 
~3~ 



m 



(m + 2Y 



+ 



m 
16 



2 (m + 1) 2 m+1 



T + 1 ) 

- 1 mj 
16 + 4 



m 



(4.16) 



(4.17) 
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and 

f(m)=4(^^ + -L-) , (4.18) 
\m + 2 m + 1/ 

where 

n i 

%) = E«- (4-19) 

™ / „ \ c_i y+i i i i 

S(n) = E (J J = 2 S ^ n - X ) + n Sl{n) + 2^ " 1} " (420) 

The same points we noted below Eq. (|4.13|) apply here and are discussed below. In 
addition, we note that the coefficient of the double pole is the same as that in Eq. ( |4.13j ), i.e. 
it is universal, and that the coefficient of the single pole contains the anomalous dimension 
of the gluon operator in addition to the terms in Eq. (|4.13|) . 



V. DISCUSSION AND CONCLUSIONS 

We now return to the three points of contention listed in the introduction. There are two 
reasons why the matrix elements of the GV operators used by Hamberg and van Neerven 
H did not vanish. One is that their operators are not BRST exact. Therefore, Theorem 2 
does not apply. One can see this by simply checking that at Born level their O a does not 
vanish on-shell when Q ^ 0. Even if they had used an operator that was BRST exact they 
would have found it had a non- vanishing physical matrix element because Theorem 2 breaks 
down when the operator is inserted at zero momentum as was shown in detail by Collins 
and Scalise [TO]. The fact that Theorem 2 breaks down at one loop for a zero momentum 
insertion means that at two loops one will see a non-triangular mixing matrix, hence the 
apparent demise of Theorem 5. By inserting BRST exact operators at non- zero momentum 
we have shown Theorem 2 works at one loop and we therefore have no reason to doubt 
Theorem 5. 

The gauge parameter dependence of the physical matrix element of the gauge invariant 
gluon operator was not investigated by the first group attempting a two loop calculation in 
the covariant gauge ||. They worked in the Feynman gauge. Hamberg and van Neerven 
noticed this gauge dependence as they calculated the 0(a s ) terms in the general covariant 
gauge. They claimed that this implied a renormalization of the gauge parameter in the 
two loop calculation to keep the coefficients of the single pole terms gauge independent. 
They then proceeded with the two loop calculation in the Feynman gauge. Hence, they 
can not make a statement on the gauge parameter dependence of the finite terms at two 
loops. A finite gauge parameter dependence at two loops can lead to 1/e gauge parameter 
dependence at three loops. This would imply a gauge dependent anomalous dimension. As 
we demonstrated, our result is independent of the gauge fixing parameter at one loop. It is 
now safe to proceed to two loops. 

In our method we put the momentum squared of the external legs to zero at the very 
beginning. This is possible because we work in n > 4 dimensions and only after extracting 
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the IR safe anomalous dimension do we let n — > 4. Because we worked with physical legs 
the LSZ residue is unity. We also notice a double and single Sudakov pole structure arising 
because we set the momentum squared of the external legs to zero in the very beginning. We 
conjecture that these factors are universal to all three point functions involving two external 
gluon legs. An all orders proof of the universality may be possible and a two loop calculation 
of the energy-momentum tensor is technically possible and would serve as a check on these 
ideas. A detailed study of the universality of double and single pole structure will also yield 
information on the finite form factor seen in Eq. (|4.13|) and Eq. ( |4.15|) . As can be seen 



from Eq. (|4.13|) it is safe to let Q — > as e > (only IR divergences remain) leaving only 



the Born level result. The same is true for Eq. Q4.15| ) after the UV counter term has been 
added. 

The fact that it took so long to verify the general theorems was anticipated in |J where, 
following the proofs of the general theorems, there are statements like 

. . . matrix elements of N(x) vanish if the on-shell limit can be taken. . . . but for 
pure Yang-Mills theory, the on-shell limit cannot be taken because of infrared 
divergences . . . 

It is only by seeing that these infrared divergences are universal if regulated dimensionally 
that we are able to overcome these problems and verify the general theorems. Our method of 
calculation has actually been applied to the two loop anomalous dimension of the nonsinglet 



quark operator |22j where the complications of GI and GV operator mixing do not enter 
and where the ct^CV and a^n/CV terms were found to be in agreement with the Sudakov 



factorization theorems 11 



To conclude, we have presented new one loop calculations that vindicate the general the- 
orems of Joglekar and Lee which were recently called into question by the work of Hamberg 
and van Neerven. In the process we have found a new method of calculating the anoma- 
lous dimension of the gluon operator. We did this by considering physical matrix elements 
with insertions of the appropriate operator at non-zero momentum with the resulting IR 
singularities being regulated dimensionally. We found that the physical matrix element of 
either the classical energy-momentum tensor or the gauge invariant gluon operator is in- 
dependent of the gauge fixing parameter. We also found that a physical matrix element 
of the BRST exact gauge variant operator that enters the full energy-momentum tensor 
vanishes. A universal Sudakov factor appears in both the physical matrix elements of the 
energy-momentum tensor and of the gluon operator. The universality of this factor and the 
UV finiteness of the energy-momentum tensor then provide another method for determining 
the gluon anomalous dimension which we conjecture applies to higher loops as well. 
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APPENDIX A: FEYNMAN RULES 

The operators Of^, O^ 1 and 0^ have the following Feynman rules as shown in Figures 
3 and 4. All momenta flow into the vertices. 

T;i;i{k x , fa, M) = -U a ^ [(A • fa) M ' 2 + (A ■ fa) M ~ 2 ] (Al) 

x [A m A M2 fa- fa- A^k lfl2 A -fa 
- A^fa^A-fa + A-faA-fag^) , 
U Zl%( k ^ fa) = -igf a ^ (A2) 

+ A /i2 (A;2 WJ A w - A^fa^) + A.fe 2 (s Mm A M - ^ 2M3 A Ml ) 
+ A M (fc 3w A M2 - A^fa^) + A-faig^A^ - ^ sm A M2 )] , 

^(Aji, *&) = -^ aiaa [A;f A^ 2 A ■ fcj + fcf A«A ■ fc 2 ] , (A3) 
W° ia2 (fci,fc2) = <^ aia2 A- fa A- fa, (A4) 

^;r 2a3 (^) = -^r ia2as a. a* a M3 . (as) 

T^ 2 2 (ki, fa, M = 2) is the Born level energy-momentum operator vertex and V^*(fa,k 2 ) 
vanishes when contracted with polarization vectors satisfying e(k) ■ k — 0. 



APPENDIX B: INTEGRALS 

In this appendix we briefly discuss some technicalities of the calculation. For the case 
of Q = all integrals are expressible in terms of Eq. ( P2|) below. The case of Q 7^ is 
more complicated and we discuss it now. After partial fractioning the amplitudes we find 
everything is expressible in terms of the integrals 
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J(a,b, c, M){ tltttMl/ y 



d n k {l,k^k^k u }(A-k) M 
(2ir) n (k 2 ) a [(k + Pl ) 2 } b [(k + p 2 ) 2 ] c 



(Bl) 



We write these integrals as linearly independent tensor structures multiplied by undeter- 
mined coefficients. We then contract both sides with the various independent tensor struc- 
tures and get a set of linear equations. These we invert and solve for the previously unde- 
termined coefficients. At this point we have expressions for the tensor integrals above in 
terms of the scalar integral J(a, b, c, M) and the tensor integrals with only two denominators, 
namely 



I(a,b, M) { ^ v} 



d n k {1, fc M , k fl k u }(A ■ k) 



M 



(B2) 



(27r) n (k 2 ) a [(k + p) 2 } b ' 
which we again decompose as above. The remaining integrals are expressible in terms of 

/(a, b, M) = tS n (-l) n/2 (p 2 ) n/2 - a - b (-l) M '(A ■ p) M (B3) 
r(a + b - n/2) T(n/2 - a + M)V{n/2 - b) 



x 



T(a)T(b) 



r(n + M — a — b) 



and, for the p\ =p\ = m 2 case, the specific cases J(l, 1,1,0), J(l, 1,2,0) = J(l,2, 1,0), 
J(2,2,l,0) = J(2, 1,2,0) and J(2, 1,1,0). In the above formulae 5 n = 7r n / 2 /(27r) n . The 
scalar J integrals listed above can be found in the literature [23] . We list the cases of interest 
here for completeness. The integral J(l, 1, 1, 0) is finite in four dimensions for m 2 ^ 0: 



J(l, 1,1,0) 



16vr 2 Q 2 



where 



$(x) 



1 



X(x) 
2Li 2 



|21n 


x — X(x) 


In 




2 





x — \(x) 


-2Li 2 


2 





x§(x) 



2 — x — X(x) 



2 - x - X(x) 
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+ T 



where Li 2 is the dilogarithm function and 



\(x) = X\Jl — A/x 

x = Q 2 /m 2 . 

From the above and properties of the dilogarithm function one can work out that 



Re[x$>(x)] 



In (x) as x — >• oo 
ln 2 (— x) as x — > — oo 



(B4) 



(B5) 



(B6) 
(B7) 



:b8) 



The other J integrals have poles in n — A. They are, dropping terms 0(n — 4), 
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J(2,l,l,0)=iS n 
J(l,2,l,0) = iS n 
J(2, 1, 2, 0) = iS n 



n 




4 




2 




n 




4 




2 





-In 
+ In 



n 



4 



(Q 2 



,2\ 



+ 2 In 

+ 7 



-m 



+ 7 



rrr 



m 



+ (m 2 - Q 2 ) In 



. » 2 



Q 2 m 2 
Q 2 + 2Q 2 \n 
1 



-m 



(B9) 
(BIO) 
(Bll) 



7 



g 2 m 6 



where p 2 = p\ = m 2 , Q = pi — p 2 , 7 is the Euler constant and the renormalization scale, 
which comes from the coupling constant implied on both sides of the formulae, is chosen to 
be fi 2 . Note that the terms in ln(/i 2 ) cancel in the final answers given in Section IV. 
Whereas for the p\ = = p\ case we need 



J(a, b, c, M) = i(-l) n/2 S n (Q 2 ) n/2 ~ a ~ b ~ c (-A ■ p 2 ) M 

T(a + b + c- n/2)r(n/2 - a - c) T(n/2 + M - a - b) 



(B12) 



r(6)r( c ) 



r(n + M - a - b 



xF 2>1 i-M,--a-c;l---M + a + b; 

\ 2 2 A-p 2 J 

where Q 2 = —2p\ -p 2 7^ 0, and F 2t x(a, (3; 7; z) is the hypergeometric function of one variable. 
Both scalar integrals Eq. ( P^ ) and Eq. ( |B12| ) can be derived by repeated differentiation of 
the M = integral or by the method of f23| . We have checked that both methods agree 
with each other and, for M = 2,4 and 6, with the standard Passarino-Veltman reduction 
techniques P5| . As stated in Section II, we use the condition A • p\ = A ■ p 2 . Thus, the 
hypergeometric function in Eq. (|B12|) reduces to products of Gamma functions. 
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